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A variant of evaluation of the theoretical values of the surface energy of quartz glass and various forms of 
crystalline silica based on analysins of the structures of these materials and experimental detennination of the 
zeta-potential of P-quartz and amorphous silica is considered. 


The surface energy of solids determines numerous pro¬ 
perties of powders: the rate of dissolution in liquid media and 
melts, the rate of chemical reactions in heterogeneous sys¬ 
tems, the sintering rate, the adsorption capacity of the pow¬ 
ders (including hygroscopic adsorption), their tendency to 
cake, the specific heat of wetting, etc. Knowledge of the spe¬ 
cific surface energy ct makes it possible to evaluate the ener¬ 
gy spent on formation of interphase surfaces in milling. 
Therefore, in solving various theoretical and applied prob¬ 
lems, much attention is given to aspects of experimental de¬ 
termination and calculation of a for various materials. 

There are several methods for determination of c for so¬ 
lids, including crystalline and vitreous silica [1]. 

A drawback of methods based on determining the sorp¬ 
tion capacity or the heat of wetting of the surface of a crystal¬ 
line body by a liquid consists in lack of reliability of the 
juvenile purity of freshly fractured surfaces, accurate know¬ 
ledge of the fracture surface relief, and accurate values of the 
interphase surface area. The method for evaluating c for 
quartz glass by extrapolation of surface tension values for 
Si0 2 melt at high temperatures (1300°C or more) and down 
to room temperature is not correct, since the surface tension 
of a melt is equal to the work done to form a unit of the melt 
surface by deformation without disturbing the liquid-me¬ 
dium continuity. The surface energy of a solid is equal to the 
work done to form a unit of its surface in destruction (distur¬ 
bance of the medium continuity), i.e., the notions of surface 
tension and specific surface energy are not equivalent. 

Methods for determining the parameter ct for amorphous 
(vitreous) silica based on surface tension data for glass melts 
extrapolated to the content of R-,0 and RO in them equal to 
zero are even less correct. 

R. Iler in his book [1] indicates a for quartz equal to 
980 x 10 _3 J/m 2 (980 erg/cm 2 ). 
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It should be noted that experimentally obtained values of 
ct for quartz and quartz glass are affected by internal defects 
and impurities existing in the structure of the solids. Thus, 
the impurity content in very pure piezoquartz and optical 
quartz glass can reach 10^ 4 — 10 5 %, or 10 _6 — 10 7 in 
fractions. In this case, the impurity content N A is expressed 
by the enormous quantity 6.02 x 10 23 x 10 ~ 7 = 6.02 x 10 16 
impurity particles per mole of silica (60 g). 

At the same time, it is of interest to analyze and compare 
calculated (theoretical) values of a for different silica modi¬ 
fications in the ideal variant (without impurities), which 
makes it possible to carry out a comparative evaluation of the 
surface properties of crystalline and vitreous Si0 2 . 

In the principal forms of SiO.,, each silicon atom is sur¬ 
rounded by four oxygen atoms, and each oxygen atom is sur¬ 
rounded by two silicon atoms. According to Poling, the de¬ 
gree of covalence of the Si - O bond constitutes 65% (ac¬ 
cording to other sources, it even reaches 80 - 89%). The 
length of the Si - O bond in different silica modifications 
varies from 0.158 (high-temperature cristobalite) to 0.166 nm 
or more, its mean value being 0.162 nm [2]. This value is 
taken as the basic value. In this case, the distance between 
the vertices of the silicon-oxygen tetrahedron (i.e., the dis¬ 
tance between the oxygen atoms in it) constitutes 0.265 nm. 
The molar-energy values for the single bond E Si _ 0 vary as 
well for different modifications of Si0 2 : 106, 102.8 (a-cri- 
stoballite), and 102.25 (quartz glass) kcal/mole [2, 3], 2 We 
take i? si _ 0 = 106 kcal/mole = 444 kJ/mole. 

Consequently, the single-bond energy in Si0 2 is equal to: 

£ si -o = £'si-o/^ = 444/(6.02 x 10 23 ) kJ/bond. 

When a solid is broken along some plane ( hkl ), two frac¬ 
ture surfaces with a certain number of lattice points that have 
broken bonds with the part of the crystal fractured off arise 

2 1 kcal/mole = 4.1840 kJ/mole. 
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(their reticular density is d™}.). If each of these lattice points 

has a certain number of broken (uncompensated) bonds 
(n br ), it follows from the law of energy conservation that the 
specific surface energy of the solid is equal to 

a hki = ~ £ si - o d-m n br - ( 1 ) 

In the case of Si0 2 , it is easiest to evaluate and n br 

for the a-cristoballite lattice, belonging to the cubic system 
(Fig. 1). The cell parameter a for this lattice is equal to 
0.716 nm. 

Theoretically, it is possible to compare the variants of 
a-cristoballite crystal fracture along the planar networks 
{001} and {111}, although it is known that cleavage in 
a-cristoballite (as in other forms of silica) is lacking or ex¬ 
tremely unclear [4]. The fractures of all silica modifications 
usually have a conchoidal structure, but this does not prevent 
certain theoretical reasoning. 

Assuming that a-cristoballite breaks along the planar 
network {001}, one of the fracture surfaces will have only Si 
atoms with two broken bonds (known as a siliconized sur¬ 
face), while the other surface will contain twice as many O 
atoms with one broken bond (an oxylated surface). The sur¬ 
face area of an a-cristoballite cell on the network {001} is 
a 2 = 0.716 2 = 0.5127 nm 2 . Statistically, there are two silicon 
atoms on this area (see Fig. 1) (one atom at the center of a 
cell face, and also 4 x 1/4, i.e., one, lattice point from the 
4-vertex lattice points). The number of silicon atoms calcu¬ 
lated per nm 2 of a-cristoballite fracture surface on the net¬ 
work {001} is equal to 2 x 1/0.5127 = 3.9, and the number 
of broken bonds is equal to 3.9 x2 = 7.8. Then, 

1 , 

CT ooi = 2 E si-o x 3 - 9 x 2 = 3.9e si _ 0 per nm . (2) 

The quantity a 001 for oxylated surfaces will be the same, 
assuming that the energy of the single bond s Si _ 0 is equally 
distributed between the two surfaces of a-cristoballite frac¬ 
ture. 

Figure 2 shows the structure of a-cristoballite projected 
on the plane {111}. The structure can be represented as a 
three-dimensional skeletal chain of silicon-oxygen tetrahedra 
[Si0 4 ] that are helically joined to each other. These tetrahe¬ 
dra in the lattice structure are located at different levels. The 
vertices of tetrahedra [Si0 4 ] of the first level 1 are directed 
upward: these are broken Si - O bonds. To preserve the elec¬ 
trical neutrality of the fracture surface, these vertices are al¬ 
ternately electropositive (=Si - 0 + ) and electronegative 
(=Si-0“). The vertices of tetrahedra [Si0 4 ] of type 2 lo¬ 
cated below the indicated level are directed downward, and 
these tetrahedra are linked to tetrahedra of a lower level 3, 
etc. (3 linked to 4, 4 to 5, 5 to 6). It can be seen that in the 
case of a perfect fracture of a-cristoballite along the plane 


001 



Fig. 1 . Elementary cell of a-cristoballite. 



Fig. 2. Projection of the a-cristoballite structure on the plane {111}. 

{Ill}, broken bonds will be observed only in tetrahedra of 
type 1 (one bond per tetrahedron), i.e., not all points of the 
Si0 2 lattice located on the fracture surface will have uncom¬ 
pensated bonds (in the perfect case). 

It is possible to distinguish an elementary loop of area 
/ x h on the fracture surface {111}, and it is easy to calculate 
that there is one broken Si - O bond per elementary loop 
(half a bond each from two tetrahedra of type 1 ). These 
loops completely fill the fracture surface and form a pattern 

similar to brickwork, with a loop shift of ^ /. The quantity 1 is 

equal to the sum of two edges of a tetrahedron (2 x 0.265 = 
0.53 nm). The loop height h — 2 x 0.265 sin 60° = 0.458 nm, 
and its surface area / x h = 0.530 x 0.458 = 0.2438 nm 2 . The 
number of lattice points with uncompensated bonds calcu¬ 
lated per nm 2 of the network {111} will be equal to 
1 x 1/0.2438 = 4.1, and 

1 , 

CT m = 2 e Si-o x 4.1= 2.05s si _o per nm . (3) 

A comparison of expressions (2) and (3) shows that the 
energy spent on breaking a-cristoballite along the network 
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Fig. 3. Model of a hydroxylated surface of silica: 1 ) micelle nu¬ 
cleus; 2 ) surface film of silicic acid; 3 ) potential-determining ions; 
4 and 5 ) adsorption and diffusion counterion layers, respectively. 


{001} is almost double the energy spent on fracture along the 
network {111}. This means that the second variant of frac¬ 
ture is highly preferable. Obviously, the presence of real de¬ 
fects and impurities in cristobalite (and in P-quartz and 
quartz glass) dramatically alters the character of the fracture 
surface (which is conchoidal rather than flat); however, the 
tendency to have a certain portion of the tetrahedra [Si0 4 ] 
with one uncompensated bond on this surface is retained. 

The following argument can be offered to support this 
point of view. When Si0 2 breaks along the plane {111}, the 
fracture surface is electrically neutral. In the case of a break 
along the plane {001}, the oxylated surface will be 
electronegative and the siliconized surface will be 
electropositive. Since the electronegativity of Si and O is, re¬ 
spectively, 1.9 and 3.6 [5], the specific surface energy of the 
oxylated fracture surface is actually twice that of the 

siliconized surface CTgij, i.e., one surface receives 2/3e Si _ 0 
from each broken Si - O bond, and the other surface gets 
l/3s Si _ 0 , and not l/2s Si _ 0 to each surface, as occurs in the 
simplified variant. A substantial difference in the specific 
surface energies Act 001 is observed between the indicated sur¬ 
faces (an affinity of the fracture surfaces), a consequence of 
which is strong cohesion and restoration of the lattice. 

Considering that s Si _ 0 = 444/(6.02 x 10 23 ) kJ/bond, 
1 kJ = 10 3 J, 1 J = 10 7 erg, and 1 cm 2 = 10 14 nm 2 , we obtain 
for a-cristoballite: 


444x 10 3 x 10 7 x 10 14 


2x6.02x10 


23 


x 7.8 = 


2875.32 erg/cm 2 * 2.875 J/m 2 ; 


= 444 xlO 3 xl0 7 xl0 14 
111 2x6.02xl0 23 

1511.38 erg/cm 2 » 1.511 J/m 2 . 


It is of interest to note that the calculated values of ct for 
the diamond lattice are ct u1 = 5620 erg/cm 2 and g 001 = 
3820 erg/cm 2 [6], It can be seen that for these materials ct 
has a certain correlation with the Mohs hardness (10 for dia¬ 


mond and 7 for SiO,). It is known that the hardness as well is 
determined by the quantity s. 

Knowing the density p of a- and P-cristoballite (2.20 and 
2.32 g/cm 3 ), P-quartz (2.65 g/cm 3 ), and quartz glass 
(2.20 g/cm 3 ), it is easy to evaluate the number of mechani¬ 
cally broken bonds per unit surface area (1 nm 2 ) for these 
materials: 

n =_ n a-ci _ = __- 425' 

P ~ Cr (Pa-cr/Pp-cr) 2/3 (220 / 232) 2 ' 3 

n =__ = _«5_ = 463 . 

P ~ q (Pp-cr/Pp- q ) 2/3 (232/2.65) 2 / 3 


„ =_12=®_ = __= 413 

q ' S (Pp-cr / Pq.g) 2/3 (232 /220) 2 / 3 

By substituting the obtained values of n into the formula 
ct = l/2s Si 0 «, we get 

CTp. cr = 1566.68 erg/cm 2 w 1.567 J/m 2 ; 

CTp. q = 1706.75 erg/cm 2 « 1.707 J/m 2 ; 

CT q g = 1522.45 erg/cm 2 * 1.522 J/m 2 . 

Apparently, the coincidence of the calculated values of ct 
for a-cristoballite and quartz glass (1.511 and 1.522 J/m 2 , 
respectively) should be indubitable, since their densities are 
equal. However, the actual situation is more complicated. 

Determination of the electrokinetic potential £, of fine 
(1-30 pm) particles of pure P-quartz and quartz glass held 
in doubly distilled water for 4-5 days, using the micro¬ 
electrophoresis method, showed that the ^-potential of parti¬ 
cles of hydrated P-quartz (- 42 mV) is lower than that of par¬ 
ticles of hydrated quartz glass (-45 mV), i.e., the actual 
electrical-surface activity of quartz glass in water is 1.07-fold 
higher than that of P-quartz. 

At the same time, in hydroxylation of surface lattice 
points of these materials, the active centers are primarily 
sites with broken (uncompensated) Si - O chemical bonds 
(i.e., =Si - O 4 and sSi - 0~). Presumably the surface of hy¬ 
drated silica particles (micelles) to a great extent inherits the 
type of fracture surface observed in pure silica modifications. 
Silicon-oxygen tetrahedra with their vertices alternately di¬ 
rected upward (toward the dispersion medium) and down¬ 
ward (toward a micelle nucleus) are located on the hydrated 
fracture surface (Fig. 3). Due to reaction with water mole¬ 
cules vertices oriented toward the liquid transform into OH 
groups, which arise according to the following scheme: 

=Si 0 + + HOH -> =Si - OH + H + ; 

=Si - 0“ + H + h> =Si - OH. 

Two other vertices of [Si0 4 ] tetrahedra are directed to¬ 
ward the viewer and into the plane of the figure, ensuring the 
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presence of a continuous surface film of silicic acid of com¬ 
position H-,Si 4 0 9 enveloping the micelle. Partial dissociation 
of this acid leads to the appearance of a negative surface 
charge on the micelles (potential-forming O particles) and a 


^-potential at the gliding interface. 

The hydrogen particles contained in nondissociated OH 
groups form an adsorption counterion layer, and the protons 
H + arising in the dissociation of OH groups form a diffusion 
layer. 

The higher ^-potential of hydrated quartz glass particles 
(- 45 mV) compared to that of (3-quartz particles (- 42 mV) 
is the evidence of a larger number of OH groups per unit sur¬ 
face area of the micelles in glass (and, consequently, a larger 
number of ruptured Si - O - Si bonds per unit surface area 
of quartz glass) compared to the surface of (3-quartz, i.e., 


actual 

n 

q-g 


> n \ 


i-q ' 


It is obvious that « aC g Ual / H p- q 


5,. s Ai 


>-q> 


whence 


actual 

q-g 


= n 


P-q 


^ = 4.63x 

Vq 



4.95 per nm 2 . 



The difference n actual - n calc = 4.95 - 4.13 = 0.82 bond /nm 2 
q-g qg ' 

correlates with the number of thermally broken bonds per 
unit surface area of quartz glass particles. Indeed, in produc¬ 
ing Si0 2 melt at high temperatures (2000°C or more), a sub¬ 
stantial number of Si - O - Si bonds are thermally broken, 
which is accompanied by uncoupling of electrons in sp 3 va¬ 
lence orbitals of Si atoms and p orbitals of O atoms. In the 
case of rapid cooling of the melt, these orbitals do not have 
time to overlap each other, and the thermal ruptures of 
Si - O - Si bonds become “frozen.” An example of such a 
rupture in Fig. 2 can be rupture of the bond between tetrahe- 
dra of types 1 and 2 along the bridge oxygen atom (the com¬ 
mon vertex). Such breaks impart a certain looseness to the 
quartz glass structure (a kind of nonequilibrium internal de¬ 
fect), which is one of the most important reasons for its ther¬ 
modynamic instability. The same reason explains why the 
energy consumed on formation of unit interphase surface in 


quartz glass (a = 1.522 J/nm 2 ) is lower than in (3-quartz 

(Cp q = 1.707 J/nm 2 ), while the overall of surface energy 
level, on the other hand, is higher in quartz glass 
(a actual = 1.825 J/m 2 ). Quartz glass arrives at the milling 

stage with a substantial number of Si - O - Si bonds in its 
structure already thermally broken (inheritance from the 
melt), which is not the case in (3-quartz. 

The present paper introduces certain corrections to the 
values of n for p-quartz and quartz glass quoted in [7]. The 
reasoning in the latter publication was based on the number 
of broken bonds per nm 2 of cristobalite surface equal to 4.6, 
which was taken from R. Iler’s book [1], A more detailed 
analysis of the a-cristoballite lattice indicated that for this 
material n = 4.1, and n = 4.6 should refer to P-quartz. The 
value « aC g Ual =5.37 for quartz glass was given in [7], al¬ 
though this parameter most likely is 4.95, and the number of 
thermally broken bonds per nm 2 of quartz glass surface area 
is equal to 0.82 rather than 0.9. However, these discrepancies 
do not disprove the fundamental approach to the calculation 
of theoretical values of specific surface energy for crystalline 
and vitreous silica. 

Experimentally found quantities a frequently have va¬ 
lues that do not coincide with theoretical-calculation data 
and even significantly differ among themselves, due to the 
complexity of the experiment, lack of guarantee of the juve¬ 
nile purity of freshly formed interphase surfaces, the impos¬ 
sibility of total exclusion of gas adsorption on these surfaces, 
and the effect of impurities contained in the silica (even in 
very small quantities). 
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